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Abstract 

We evaluate the size of isospin breaking corrections to form factors / and g of the i^^4 decay pro- 
cess TT^ir^i^i'i which is actually measured by the extended NA48 setup at CERN. We found 
that, keeping apart the effect of Coulomb interaction, isospin breaking does not affect modules. 
This is due to the cancelation between corrections of electromagnetic origin and those generated 
by the difference between up and down quark masses. On the other hand, electromagnetism affects 
considerably phases if the infrared divergence is dropped out using a minimal subtraction scheme. 
Consequently, the greatest care must be taken in the extraction of tttt phase shifts from experiment. 
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I. INTRODUCTION 



Measuring the quark condensate remains the main concern for physicists of non pertur- 
bative quantum chromodjTiamics. The purest process allowing a direct measurement of this 
parameter is vrvr scattering. Information concerning the latter can be obtained from the 
rescattering of two pions in the final state of pionium Q], ^ Svr [2j, or K^^ decays. Let 
5l be the phase of a two-pion state of angular momentum / and isospin / and consider the 
charged KpA decay process 

K^{p) ^ ^^{pi>-{P2)t{p,)u,{p,) , (1) 

where the lepton i is either a muon /i or an electron e, and v stands for the corresponding 
neutrino. In the isospin limit, the decay amplitude A for process (0) can be parameterized 
in terms of three vectorial (F, G, and R) and one anomalous {H) form factors: 

[(pi + P2rF + (pi - p^rc + (p, + p^YR] 

Mk± 

(2) 



where Vus denotes the Cabibbo-Kobayashi-Maskawa flavor- mixing matrix element and Gp 
is the so-called Fermi coupling constant. Note that form factors are made dimensionless by 
inserting the normalizations, M~^\ and M~^^. The fact that we have used the charged kaon 
mass is a purely conventional matter and corresponds to the choice of defining the Isospin 
limit in terms of charged masses. 

Form factors are analytic functions of three independent Lorentz invariants, 

St, = {Pl+P2f, Si = {pi+Puf: (3) 

and the angle formed by p^, in the dipion rest frame, and the line of flight of the dipion 
as defined in the kaon rest frame 14,0]. In the following, we will be interested only in two 
form factors, F and G, and consider the partial wave expansion, 

F = fs{s^, si)e<^^'-^ + fp{s^, se) cos9^e<^'-^ , (4) 
G = Ms., s,)e^^'(^-) + goisn, s^) cos O^e''"^'-^ , (5) 
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where a convenient parametrization of fs, fp, gp, and qd in the experimentally relevant 
region has been proposed in Ref. 0]. 

The currently running NA48 experiment aims at measuring form factors for Ki^ decay of 
the charged kaon with an accuracy better than the one offered by previous measurement 0, 
Isj. The outgoing data on form factors contain, besides a strong interaction contribution, a 
contribution coming from the electroweak interaction. The latter breaks isospin symmetry 
and is expected to be sizable near the tttt production threshold j^. In order to extract 
TTTT scattering parameters from the NA48 measurement, the isospin breaking correction to 
form factors should therefore be under control. In this direction, we recently published 
analytic expressions for F and G form factors calculated at one-loop level in the framework 
of chiral perturbation theory based on the effective Lagrangian including mesons, photons, 



Ul to 



and leptons [l^. In the present work, we will use the method proposed in Ref. 
split analytically the isospin limit and isospin breaking part in form factors, allowing a first 
evaluation of isospin breaking effects in charged decays. 

II. A BRIEF REVIEW OF THE METHOD 

We shall start things off by the general form of the decay amplitude for process (P) as 
dictated by Lorentz covariance, 

GfV* 

^ [{Pi + P2 ff + {Pi - P2 f9 + {Pi + PuTr] 7^ 



Mk± 

The quantities /, g, r, and h, will be called the corrected Kg^^ form factors since their isospin 
limits are nothing else than the form factors, F, G, i?, and if, respectively. The tensorial 
form factor T is purely isospin breaking and has been calculated at leading chiral order in 
Ref. |10|]. The corrected form factors as well as the tensorial one are analytic functions of 
five independent Lorentz invariants, 5,^, S£, 6^^, 6*^, and 0. 6^ is the angle formed by p^, in the 
dilepton rest frame, and the line of flight of the dilepton as defined in the kaon rest frame. 
is the angle between the normals to the planes defined in the kaon rest frame by the pion 



pair and the lepton pair, respectively. Let us denote by SF and SG the next-to-leading order 
corrections to the F and G form factors, respectively. 



V2F, 



9 = ^^{l + SG]. 



V2Fo 

The analytic expressions for 6F and 6G were given in We shall distinguish between 
photonic and non photonic contributions to 6F and 6G. The photonic contribution comes 
from those Feynman diagrams with a virtual photon exchanged between two meson legs or 
one meson leg and a pure strong vertex. Obviously, this contribution is proportional to e^, 
where e is the electric charge, and depends in general on the five independent kinematical 
variables, s^r, se, Ott, 6e, and (p through Lorentz invariants like {p2 + Pe)"^, say. The non 
photonic contribution comes from diagrams having similar topology as the ones in the pure 
strong theory with Isospin breaking allowed in propagators and vertices. This contribution 
generates Isospin breaking terms proportional to the rate of SU{2) to SU{3) breaking. 



m 



rriu + rud) , (6) 



4 rUs — m 2 
and to mass square difference between charged and neutral mesons, 

= M^2± - M^o = 2Zoe^F^ + 0{p^) , (7) 
Ak = Ml± - Ml, = 2Zoe^F^ - Bo{md - + 0{p') , (8) 

or equivalently, (m^ — mu)/{ms — m), Zqc^, and — m^. The kinematical dependence is 
on three Lorentz invariants, {pi +P2Y1 {p~Pif'i and {p — P2Y which represent respectively 
the dipion mass square, the exchange energy between the kaon and the neutral pion, and 
that between the kaon and the charged pion. In terms of independent kinematical variables, 
the preceding scalars are functions of s-„, si, and cos6'^. 
It has been noted in Ref. that for 

Si = m\ (9) 

the photonic contribution neither depends on Qi nor on and, consequently, it can be written 

as 

photonic contribution = <^(sjr) + "(9(5^) cos 6'^ , (10) 



where (; and are analytic functions of s-j^. Note that, to the order we are working, that is, 
to leading order in isospin breaking, the power counting scheme we use dictates the following 
on-shell conditions to be used in the argument of <; and 



p-^ ^ Ml ^ Boirris + m) , pj ^ pi ^ ^ 2Bom . 



(11) 



With respect to the nonphotonic contribution, it depends on St^, {p — Pif, {p — P2Y and 
masses through one- and two-point functions. In order to split strong and electromagnetic 
interactions in one-point functions we use the formula 

1 1 



167r2 Ml 



A{Ml) 



(12) 



where P denotes a pion, tt, or a kaon, and Ap the difference. 



Ap = Ml± - Mlo 



(13) 



Concerning the splitting in two-point functions B{pi,mo,mi), we have to expand exchange 
energies in powers of the fine structure constant a and md — rriu- We then inject the obtained 
expansion in the expression of B{pl + S, ml -\- 5q, m\ -\- 5i) where 5, 5q, and 5i, are leading 
order in isospin breaking, 

5 , 5o, 5i,^ 0{a,md-mu) . (14) 
The final step consists on expanding two-point functions to first order in 5, and 5i, 



B{pl + 5,ml + 5,,m\ + 5,) = B{plmlml 

1 , I ''^0 \ /2 2 2\/2 2 2\ 

hJ 



327r2p2 
1 

■32^ 
1 



In 



rrir 



So 
Si 



2p\ + (ml — ml) In 



mt 



mi 



327r2pf 

+ [{ml - mlf - pl{ml + ml] t{pI, ml ml) ^5, 



with T a generic integral defined by. 



Jo ^"^0 + (1 ~ x)mf — x{l — x)p\ 



(15) 



(16) 
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Putting all this together, form factors for Kg^ decay of the charged kaon can be writ- 
ten in the following compact form which shows explicitly the splitting between strong and 
electromagnetic interactions, 

X {Sn, {P - (P - P2)^ {P2 + Ptf, ■■) = 

Mk± 



\/2Fo 



[l + U'^{s^) + V^{s^)cose^] , X = f, g, (17) 



where, 



+W^2e^ + ly; ^ , W = U ,V , (18) 



V3 

are analytic functions of s^^. If one makes the following substitutions. 



2Z^e'Fi , (19) 
A;, 2Zoe'F^ -^(Ml- , (20) 



then, equations (fTTjl and (fTH|) read. 



e 



W"" = W + We^ + W' — , (21) 
W: = + 2Z^F^{W: + Py^) , (22) 
= W: - 4{Ml - Ml)W-^ . (23) 

The aim of the present work is to determine the U functions corresponding to / and g form 
factors for Ki4^ decay of the charged kaon. 

III. ISOSPIN LIMIT 

We have 

-QU + f {2Ml^ + Ml± + M^) (m2± - M|.±; 

+ [I6(s^ - 2M2±)Li + 4(M|± - m] + s^)L2 
^0 

+ (M|.± - SM^i - + 5s,) L3 - 2(2M|.± - 7M^i)U + m]L^] 
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8F2 



5 - ^ [Ml. - 2Ml.) + i (M,^. - Ml4' 



-- (2Mi, - Ml^ - Ml) {MU - Ml.] 



8Fi 



3 + - (2M|^± + SM^i - 3M2) 

V'jr 



tl 



[Ml^ - Ml.) {Ml - Ml.] 



ml) 



+ 



[Q{2s^-Ml.)B{s^,Ml.,Ml, 



+6M^.B{s^, Ml Ml) + 9s^B{s^, MI.,MI.)] 
+1 {M^.. - Ml.)' + ^ {M^^. - Ml.)' 



B{U,M^.,Ml.) 



8F| 



2M|^± - M^. + 3Ml - - {AMI. + Ml. - SM^) Ml. 

V-TT 



+ ^ {Ml. - Ml) Ml - I {Ml. - Ml.) {Ml - Ml.) 



[l2Ml. + 21Ml. + 3Ml 



B{U,Ml,Ml.), 



3847r2F2 



-4s, - 2i, - - (2M|^± + Ml. + Ml) {Ml. - Ml.) 
- [(^K± -mj + s^)Ls + 2{Ml. + 2Ml.)U - mjLg] 



24F2 



24F2 



tl 



{AMI. - ^Ml. - Ml) 



+ ^ {2Ml. - Ml. - Ml) {Ml. - Ml.) 



A{Ml.) 



1 + i_ {Ml - 3Ml.) + 1 {Ml. - Ml.) {Ml - Ml.) 



A{Ml) 



[2{s, - AMl.)B{s,, Ml., Ml.) + {s^ - AMl.)B{s^, Ml., Ml.)] 



AFi 



[Ml. 



Ml. - 



1 



B{t^,Ml.,Ml.) 



8F2 



(25) 



In the preceding expressions, we used the notation 



(26) 



Note that in the isospin breaking correction, the same expression holds for with the 
replacement Mp± — > Mp. 

IV. NON-PHOTONIC CORRECTION 



The correction due to e reads: 



3 + f (M^-Mi) 

t'-TT 



167r2 



+ 
+ 



-^+^Ml + ^{Ml-Ml) 



2 \2 



/l(M^) - - {Ml - Ml) {Ml - Ml)A{Ml) 



3 - ^ Ml - I (M,2 - Mi) (M,; - Af 



A{Ml) + 12(5, - M2)5(5,, Ml, Ml) 



+2 

+ 



4M25(s,, Ml Ml - 4(3s, - AMl)B{s^, Ml, M^ 

2 I r 7i/r2 oj. /'7i/r2 7i/r2\2 



2\3 



+ 5MI -3U-^ {Ml -MlY-^ {Ml - M^) 



B{U,Ml,Ml) 



-2(5Mi + - 3i,) + - {3Ml - bM^) {Ml - M^) 



+i (3M2 - Mi) (M^ -Ml)+'^ {Ml - Ml) (M^ - Ml)^ B{U, M^ M^) j(^7) 



-1 + i (M,2 - Mi) 

L-rr 



167r2 



+ 



2\2 



1 + -Mi + -(M,2-Mi) 



f f2 



A(M,2) 



1 



AMI - - M^ + - (M^ - M^) (M^ - Mi) 



^(Mi) 



+ 



-1 + 1 (2Mi - ZMl - Ml - 1 {Ml - Ml) {M^ - Ml) 



A{Ml 



2\3 



Mi - 3M2 + 2t^ + - {Ml - Ml) 



B{t^,Ml,Ml) 
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+ 



-2{3Ml + Ml- - 2t^) - - {AMI - ' ^k) 

+i {AMI + M^^) {Ml - Mi) + I {Ml - Ml) {M^ - M^f B{U, M^ M^) j(28) 
That due to M\ - M\ is equal to: 



3847r2F2 



2\2 



6 - (Mi - IMl) + ^{Ml- Ml) 



24 



24F2 



9 + ^ (13M2 -Ml)-^ {Ml - 2Ml) {Ml - Ml) 



mi) 

M2 



24F2 



3847r2F2 



9 + -(M^2-Mi) A(Mi) 
6-l(5Mi + M^) 



3 1 



8-^0 ^TT 



2\3 



+1 {Ml - 2Ml) {Ml - Mi) - 1 {Ml - Ml) 



In 



Mi 



3847r2F2 



Mi 



7M,2 + + - (7Mi - 3M2) {Ml - Ml) 



2\4 



-1 (Mi - 3M,^) {Ml - Mi)2 + UmI- Ml) 

5 



T{U,Ml,Ml) 



32n^F^ 



{s^ - Ml) t{s^, Ml Ml) + ^ B{s^, Ml Ml 



2Fi 



l2Fi L 



\B{s^,MlMl)-AB{s^,Ml,M^)\ 



3 + (lOMi - 13M,2) - |- {Ml - Miy 



B{U,MlMl) 



l--^(M^^-Mi) 



B{t^,MlMl), 



3847r2F2 



2\2 



10 - - (Mi + bMl) - - {Ml - Ml) 



tl 



2AF^ 



22 2 
U--Ml + -{Ml- 2Ml) {Ml - Ml) 



1 1 



24F2 



7+ {Ml -Ml) 



mi) 

Ml 



ml) + J| ^ A{M'^) 



7 (7Mi + 3M^) 



SSAti^F^ 

+i (Mi - 3M3 {Ml - Mi) - 1 {Ml - Ml) 



2\3 



SSAti^F^ 



"TT 

r2 / 71 ^2 



In 



Ml 
Ml 



IQMl -7t^-- Ml{Ml - Ml) 

L'Tr 



(29) 



9 



''IT ''n 



8Fi 



23 + i {IIM^ - SMl) + I {M^ - MlY 



B{U,MlMl) 



B{U,MlMl). 



(30) 



The one due to — M^o is given by: 



^ 3847r2F2 



18 (20M|. - IIM^ - 9M^) 



6 (22Af|. - 7Ml - 9M^) 

t-jr 



- {%Ml - Ml - ZMl) {Ml - Ml) 



mi) 

Ml 



1 1 



4F2 U 



3847r2F2 



1 + ^ (M,2 - Ml) 
6-l(5Mi + M^) 



ml) 



1 



+^ (Mi - 2Ml) {Ml -Ml)-^ {Ml - Mlf 



''TT 



In 



Ml 
Ml 



1 



2Ml -Ml + 3M^ + - (5M2 - - 4Mi) M 

In 



+r - K) K - I {^1 - Ml) {Ml - Ml 



2 \2 



1 



TT 

s^t{s^, Ml, Ml) 



Ml 
Ml 



llMl - 5Ml -6U + - {3Ml + 5Ml) {Ml - Ml) 



3847r2F2 

-i {3Ml - 5Ml) {Ml - Ml)' + 1 {Ml - Ml)' 

''IT ''TT 

2Ml - Ml + 3M^ - ^ {Ml -MlY-^ {QMI - Ml - bM^) Ml 



2567r2F2 



r{t^,Ml,Ml) 



+^ {2Ml + Ml - 3Ml) - ^ {Ml - Ml) {M^ - Ml)' 



tl 



ti 



12Fi 



8F^ 



2\2 



3 + - {2Ml - 5Ml) - - {Ml - Ml) 



tl 



t{u,mImI) 

B{U,Ml,Ml) 



4 - i {6Ml + Ml - + 1 {Ml - Ml) {M^ - Ml)] B{U, M^ Mj^X) 
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3847r2F2 



14 - - {AMI - - M^) 



+ 



2 1^4 



1 1 



1 - - {Ml - Ml) 



1 1 



A{Ml) + ^-{Ml- Ml) A{Ml) 



24F2 



10 - - (6M| - bMl - 3M^) 



-- {ml - Ml - ml) {Ml - Ml) 



7 - - {7Ml + 3Ml) 



A{M\ 



Ml 



3847r2F2 



2\3 



+ ^ {Ml - 3Ml) {Ml -Ml)-- {Ml - Ml) 



''TT 



In 



Ml 

Ml 



2567r2F2 



2 - ;^ {2Ml + Ml + Ml) - ^ {2Ml + Ml - SM^) M; 



1 (2M|, - 3M^ + Ml) Ml - I {Ml - Ml) {Ml - Ml 



2 \2 



tl 



M2 



3847r2F2 



14Mi - 4M2 + ^ {Ml + Ml) {Ml - Ml) 



2 \4 



''TT 



(Mi - 2Ml) {Ml - Ml)' + - {Ml - Ml) 



T{t^,Ml,Ml) 



2567r2F2 



4Mi + Ml - Ml -2t^-- {Ml - Ml) Ml 

Lit 



-1 {Ml - Ml) {Ml - Ml)Ml + 1 (4Mi - Ml - 3Ml) {Ml - MI)MI 



tl 



+^ (Ml - Ml) {Ml - Mlf 



-{t^,Ml,Ml) 



B{s^,Ml,Ml) 



12Fi 



2 \2 



5--Ml--{Ml-Ml) 



tl 



B{t^,Ml,Ml) 



SFq t-^ 



4 



2Mi + MI-MI-- {Ml - Ml) {Ml - Ml) 

V-jr 



B{t^,MlMl). 



V. PHOTONIC CORRECTION 



We have 



327r2 
1 



9 + 21n 



m 



m 



+ 41n 



Ml 



+ 2 In 



Ml 



18 



{2AKi - 2QAK2 - 16/^5 - 88/^6 - 36/s:i2 + 120Xi + OXg) 



11 



1 A{mj) ^ , A{M'^) ^ A{Ml) 



ml, 



Ml 
1 



Ml 



[A{Ml) - A{Ml)] 



2 X{t^,mlMl) 
+ ^ (Mi + 2M^) +Ml- l2Ml - bml + 2U 



mj Ml - 2Ml 
+ 2 A(t,,m2,M2) 



- 2t^ 



1 + 



2U 

{Ml + 5M^ - ?,mf) 

B{0,mlMl) 
B{mlO,Ml) 



+ 
+ 



2 + 



A(t^,m2,M2) 
{Ml - Smj + 3s^) 



ml{Ml-u] 



me 



B{mj,Q,mj) 



+ 2 



UmfM^ 



\{U,Ml,Ml) X{U,mlMl) 



ZmlMl 



[-Ml{Ml - AMI) - m^Ml + SM^^ - m|) 



X^{U,mj,Ml) 
+ {Ml - AMI - ^]%] ] B{Ml 0, Ml) 



+ 



2Ml 



\{s^,mj,Ml) 
AMI 



(M: 



2mi 



X{U,Ml,Ml) ^^'^ 



{Ml-ml-t^) 



B{Ml,0,Ml) 

^[-'^ X{s.,:iMl) [^^^(^i--) 
-m^SMi - 3m| + Ss^] } B{s^, ml Ml) - B{s^, M^ Ml) 

-l[Q+liMl-Ml-3mj) 

+ X{U,M^,Ml) [^^^ ~ 2M,')(^' - ^k) + {Ml + 2MI)U] 



mo 



[2{Ml - 13Ml - 2mj) 



X{U,mlMl) 

{Ml + 5Ml - 3m|) + ^{Ml + 2Ml) 

-M'^(Ml - 4A'Q + (Ml - ml - 5mJ)iJ } i?(i„ M|, Af 

+ A(t,.<M3 (Ml-m\ + t,)Bit„,nlMl) 
+2{Ml + mj- s^)C{ml s^, M|, m^, m|, M|) 



+2(2M2 - s^)C{Ml, s^, Ml, ml, M^ Ml) 



12 



^Jm^M^ C(Mi, s^, ml 0, M^, M^) 



X{s^,mj,Mj.) 
+m? i 1 + 



{Ml - m^U] ] C{ml 0, m], 0, ml M|) 
1 



-MlMl 



-mK2M| - m| + s^)] } C(s^, s^, 0, ml M^, M^) 



777 

+m^2(Af|, + 7 Ml + 3m^2) + _! + SM^ - Sm^ 



-2(Mi - 3M,2)i, - 
+ ^{Ml + 2Ml) 



^ {2Ml + 7Ml) 



C{U,U,0,mlMl,Ml) 



Ml 



{2{Ml - Zml){ml - Ml) 
[mj{2Ml + - ml - 2mjMlMl 



X{U,mlM^) 
3m^ 
~X{t^,mlM^) 

-MliMl + ml + ^M^{Ml - ml - mj{2Ml - AMI - mlU 

+MliMl - AMDu] } C{Ml U, ml 0, Ml Ml) , 
1 



327r2 



9 + 21n te ) + 41n ( ^ 1 + 21n 



\m; 



m: 



Ml 



Ml^l 



(33) 



{2AKi + 24^2 - 144/^3 - 72Ki + ?,2K^ - AQKq - 36/^12 - 24Xi + 9Xq) 

18 

1 Aiml ^AiMl) AiMl) 1 y,,,,., .,,.2m 



4{' 

M2 



3m? 



XiU,mlM^) [t 



!^{Ml + 5Ml-3m',) 



^ {Ml + 2Ml) - {Ml - 12Ml - 5mj + 2U) 



t. 



1 [3 

2I X{s^,mlMl^ 
^rr^{Ml - 2Ml + m| - 2t^ 



B{0,mlMl) 



+ 



2 + 



X{U,mlM^) 
2ml 
X{s^,mJ,Ml) 



{Ml + ml - s^) 
1 + 



l2mlMl 



X{U,mlM^) 



B{mlO,Ml) 



{Ml - ml + s^) 
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3m« 



+ <! -4 + 2 ^ 



SM^ - 3s^ 
4M2 - 



MlMl 



X{t^,MlMl) 



X^it^.mlMl) 
-{Ml 



[Ml{Ml - AMI) + r4{Ml + ?,Ml - m|) 



+ 



1 - 



• AMI - 
Ml 



%]]B{MlQMl) 
(2M| -ml- 2s^) 



X{s^,mj,Ml) 

(Mi + - U) 



X(U,MIMI) 
-2I1 + 



B{MlAMl) 



X{s^,mlMl) 



[mj{2Ml - mj) - M^ + (M^ + mj)s^] ] B{s^, mj, M^) 



-4 ^4M^ J ^'^) + 2X{s^,mlMl) ^^"^ 



rrio 



[-2{M 



K 



X{U,MlMl) 
IZMl - 2mf) 

Ml 



[Ml - Ml - U) 



XiU,mj,M^) 
+^ (Mi + 5Ml - 3mj) - (Mi + 2Ml) 

-3M2(Mi - 4M2) + 3(Mi - AM^ - 5mj%] } B(U, M^, Ml) 
3m? 



^ Xit^M.Ml) {^l-ml^K)B[KM.Ml) 
+2(Mi + m] - s^)C{ml s^, Ml, m% ml Ml) 
+2(2Ml s^)C{Ml s^, Ml m% M^ Ml) 



-ml { 2 + 



[m,2(Mi + 7M^-m,2) 



A(i^,m|,M2) 
+M2Mi - (Mi - m|)i,] } C{ml 0, m|, 0, m|, Mi) 



-^Cis^,s^,0,mlMl,Ml) 
+ ^{3-^(m|-M,^) + 



1 



A(t^,m2,M2) 



[-M^^M 



+m|(3Mi + 5M2 - 7ml ^ 7^ (^i + ^^l - 3m^) 
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'TT 



-4M^M| + m^^{2Ml - AMI - ^f)^- 

-Ml{Ml - AMl%] ] C{Ml ml 0, M^, M|) . 




(34) 



VI. RESULTS 

We shall proceed to the numerical evaluation of isospin breaking corrections. To this 
end, we must handle all types of singularities encountered in our expressions. These are 
of three types in general: ultraviolet, infrared, and Coulomb. Although our expressions 
are ultraviolet finite, they are infrared divergent. We showed in Ref. ^] that the latter 
singularity is canceled by the emission of a real soft photon at the level of differential decay 
rate. Since we are interested in measuring form factors, a subtraction of infrared divergence 
at this level is needed. There are infinitely many choices to do so. We shall choose the 
simplest minimal subtraction scheme consisting on dropping out, from the expression of 
form factors, Inm^ terms only. Finally, Coulomb interaction between charged particles 
induces singularities due to a photon exchange between: 

(1) the kaon and a pion. This occurs at t^, = {Mk ± M^y or 



n 



AMI - {Mk ± 2M^f + m] . 



(35) 



Hence, the singularity is situated outside the allowed kinematical region, 



AMI ^ 



s^ ^ {MK-m(,y , 



(36) 



from the left. 



(2) the two pions. This occurs at 



TT 



0, AMI 



(37) 
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The former value represents a pseudo-threshold and is situated outside the allowed 
kinematical region from the left. The latter value is a normal threshold and is situated 
at the lower bound of the allowed kinematical region. The corresponding singularity 
is of great experimental importance for the present work and we will study it further 
in the following. 

(3) the kaon and the lepton. This occurs at 

= {MK±mey. (38) 

The pseudo-threshold is situated at the upper bound of the allowed kinematical region. 
The normal threshold outside the latter from the right. 

(4) a pion and the lepton. This occurs at = (M^r ± fngY or 

= {Mk - mif + 2mi{MK T 2M^ - me) . (39) 
Hence, the singularity is situated outside the allowed kinematical region from the right. 

Let us return to the Coulomb interaction between the two pions and shift the value of 
from 4M^ by an infinitesimal positive amount 

= A{Ml + g") . (40) 

We then expand our expressions in powers of q. The Coulomb singularity shows up then as 
poles in the p-plane. In order to obtain finite (regularized) results, we simply remove these 
poles allowing the numerical evaluation of form factors. 

A. Input 

We shall use the following numerical values for the various parameters j^: 
(1) the fine structure constant, 

a = 1/137.03599976(50) , (41) 



corresponding to the classical electron charge e = yivra; 
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(2) the masses of the charged leptons, 

me = 0.510998902(21) MeV, = 105.658357(5) MeV ; 

(3) the masses of the hght mesons, 

M^± = 139.57018(35) MeV, Mk± = 493.677 ± 0.016 MeV, 
= 547.30 ±0.12 MeV, Mp = 771.1 ± 0.9 MeV ; 

(4) the quark masses and condensates, 

= 134.9766(6) MeV , Mk = 495.042 ± 0.034 MeV , 
e = (1.061 ±0.083) x 10"^ 

(5) the low-energy constants in the strong sector, 

L[ = (0.46 ± 0.24) X 10"^ , = (1.49 ± 0.23) x 10^^ , 
LI = (-3.18 ± 0.85) X 10"^ , LI = (0.53 ± 0.39) x 10"^ , 
LI = (5.5 ±0.2) X 10-^ 

(6) the low-energy constants in the electromagnetic sector, 

K{ = -6.4 X 10-3 , = -3.1 X 10-3 , = 6.4 x 10"^ , 
Kl = -6.4 x 10-3 , k; = 19.9 X 10-^ , = 8.6 x 10-^ , 

Kl^ = -9.2 X 10-3 , 

with an error of ±6.3 x 10-3 assigned to each of them; 

(7) the low-energy constants in the leptonic sector, 

I Xi 1^ 6.3 X 10-3; 

(8) the coupling of axial currents to the vacuum, 

57.40 ^ Fo ^ 67.53 ; 

(9) the charged pion decay constant and electromagnetic mass, 

= 92.419 ± 0.325 MeV , Zq = 0.805(1) . 
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FIG. 1: Radiative correction to the real part of the first term in the partial wave expansion for / 
form factor under the assumptions S£ = mf = rrig, Fq = 67.53 MeV. The plain curve represents the 
one-loop correction in the absence of isospin breaking. The dashed curve gives the isospin breaking 
correction of order 0{a, nid — niu)- The infrared divergence has been removed applying a minimal 
subtraction scheme. Error bands come exclusively from the uncertainty in the determination of 
low-energy constants and have been developed in quadrature. 

B. Form factors 

Using the preceding input parameters, we drew (Fig. ^ the curve of the variation of 
one- loop level correction to the real part for / form factor as function of s.,^. In Fig. |21 we 
drew the isospin breaking correction to the same quantity and compared the contributions 
of orders 0{a) and 0{md — mu). 

The same has been done for the real part of g form factor in Figs. 01 and 0] 
After removing infrared singularity and Coulomb poles, the NA48 experiment should 
measure what we will call subtracted form factors. The corresponding modules are found to 
be 

fsisw) = 1 + 3? [/-^(stt) + subtraction , (56) 

gpisn) = 1 + 3? f/^(s,r) + subtraction , (57) 
M 

subtraction = ^ If (s^ = 4M^) 

16 g 



87r2 



ln(m: 
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FIG. 2: Isospin breaking correction to the real part of the first term in the partial wave expansion 
for / form factor under the assumptions = = ml, Fq = 67.53 MeV. The infrared divergence 
has been removed applying a minimal subtraction scheme. Error bands come exclusively from the 
uncertainty in the determination of low-energy constants and have been developed in quadrature. 




FIG. 3: Radiative correction to the real part of the first term in the partial wave expansion for g 
form factor under the assumptions se = mj = rrig, Fq = 67.53 MeV. The plain curve represents the 
one-loop correction in the absence of isospin breaking. The dashed curve gives the isospin breaking 
correction of order 0{a, rrid — ruy). The infrared divergence has been removed applying a minimal 
subtraction scheme. Error bands come exclusively from the uncertainty in the determination of 
low-energy constants and have been developed in quadrature. 
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FIG. 4: Isospin breaking correction to the real part of the first term in the partial wave expansion 
for g form factor under the assumptions = = ml, Fo = 67.53 MeV. The infrared divergence 
has been removed applying a minimal subtraction scheme. Error bands come exclusively from the 
uncertainty in the determination of low-energy constants and have been developed in quadrature. 

where \i{ar gument) is a logical function equal to 1 if argument is true and to if argument 
is false. 



Phase shifts 



The S-wecve iso-scalar and P-wave iso-vector tttt phase shifts are given by 



7r± 



(59) 
(60) 



respectively. The imaginary part of the first term in the partial wave expansion of form 
factors reads: 



1 - 



2 ln((7^) - In ( — i 



m: 
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FIG. 5: The imaginary part (in radians) of the first term in the partial wave expansion for / form 
factor under the assumptions si = mj = ml, Fq = F^^ = 92.419 MeV. The infrared divergence 
has been removed applying a minimal subtraction scheme. The plain curve represents 5g(s^). The 
dashed one includes isospin breaking effects. 



1 + 



3a; 

+- 

a 



3a , ^ 

+-I1 + 



M2 -m? 



2M^ 
Ml 



U J - {me + M^y^t^ - (m^ - MJ^ 
1 - ln((T^) 



(61) 



2 



nit 



(62) 



K J - {me + M^Y^K - ijni - M^f 
If the infrared divergence is removed using a minimal subtraction scheme the imaginary 
part for / form factor takes the shape of Fig. El In Fig- El "we compared the size of each 
contribution to the isospin breaking part of the same quantity. Finally, the imaginary part 
for g form factor is sketched in Fig. [71 Note that the isospin breaking part is purely of order 
0{a). 

For experimental purposes, we define the subtracted phase shifts as: 



^s{s,,) 



6p{s^) 



'^Uf{s^ 
a 

~8 ~ 

a 
2 

^UHs^ 



[3Zo + 41n(2f?)]If(s^ = 4M2) 



+ 



2Ml 



In(m^) , 



(63) 
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FIG. 6: Isospin breaking correction to the imaginary part (in radians) of the first term in the partial 
wave expansion for / form factor under the assumptions = m| = m^, = -^tt = 92.419 MeV. 
The infrared divergence has been removed applying a minimal subtraction scheme. 
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FIG. 7: The imaginary part (in radians) of the first term in the partial wave expansion for g form 
factor under the assumptions = m| = rn^, Fo = F^^ = 92.419 MeV. The infrared divergence 
has been removed applying a minimal subtraction scheme. The plain curve represents SKstt). The 
dashed one includes isospin breaking effects. 

z g 

a ( ^ 2Ml\ 1 , . 2^ 



+ - 1 ^ - In(m^) . (64) 
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VII. CONCLUSION 



In this work we made the sphtting between strong and electromagnetic interactions in 
K^i decay of the charged kaon, Ti^n^i^Ui. Our expressions were evaluated at the 

production threshold for the lepton pair, S£ = mj. Thanks to this assumption, a partial 
wave expansion of form factors with exactly the same structure as in the pure strong theory 
was possible. The imaginary part of such an expansion involves the S-wave isoscalar and 
P-wave isovector vrvr phase shifts, 5o(s^) and 5{(s,r), respectively. These can be related to 
TTTT scattering lengths via Roy equations. In their turn, scattering lengths are sensitive to 
the way Chiral symmetry is spontaneously broken. Consequently, a theoretical study of the 
process in question including all possible contributions is imperative. We gave here the first 
analytic and numerical evaluation of the isospin breaking contribution. This would allow 
the extraction of 5g(s^) and SKs-,^) from the experimental measurement of form factors. Our 
results can be summarized as follows: 

• Isospin breaking affects modules of form factors only by the effect of Coulomb interac- 
tion between charged particles. The one between the two pions is of great importance 
and induces a singularity at s^, = 4M^. We gave the residue of the pole in the present 
work. 

• The effect of isospin breaking on the imaginary parts of form factors is considerable if 
the infrared divergence is removed using a minimal subtraction scheme. We gave here 
all analytical expressions for the imaginary part including the finite part, the infrared 
divergent part, the singular part with the residue of the pole. 

Our results are of great utility for the interpretation of the outgoing data from the upgraded 
NA48 experiment at CERN. 



APPENDIX A: LOOP INTEGRALS 



We use dimensional regularization and adopt the MS subtraction scheme 



A 



1 



327r2 



4 - n 



+ 1 - 7 + ln(47r) 



(Al) 



where n is space-time dimension and 7 the Euler constant. All the technical material 
necessary for the calculation of one- loop integrals is given in the appendix of Ref. 
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It is convenient to take the following notations: 



ap = Wl 



app 

\^/\U,MlMl) 

Xi 



AMI 



ap — 1 

(7P+1 ' 

y/{MK + -s^- y^{MK - - 
^J{Mk + m^)2 - + ^J{Mk - - ' 
^/U - {me + M^y -yft^-{mi- M^f 
^/U - {me + M^f + V^TT - {rm - M^f ' 
v/(M^ + M^)2^ - v/(M^ - MkY - U 
y(M, + Mk)^ -U + V(M, - Mk^ - U ' 
(m£ - AfA-)2 - s^\J{me + Ma')^ - , 

y/U - {me - M^y^yu - {me + M^f , 

y/{M, - MkY - W{M. + MkY - , 



X{U,MlMl)+4U{Ml 
X'/'{U,M^,Ml), 



m 



■e) ' 



(A2) 
(A3) 
(A4) 

(A5) 

(A6) 

(A7) 
(A8) 
(A9) 
(AlO) 
(All) 



1. A integrals 

The one-point function reads: 

A{m^) 



m 



-2A 



1 ^ / 
m 



167r2 



where /i an arbitrary scale with mass dimension. 



(A12) 



2. B integrals 

We need the following two-point functions: 
S(m|,0,m|) 

B{M^,0,M^) 



B{Ml,0,Ml) 
B{0,mlMl) 



mi 

Ml 
A{Mf 



K) 



Ml 



+ 



+ 



167r2 ' 
1 



A{m1) 



+ 



1 

16^' 
1 



Ml 



nil 



167r2 Ml 



In 



mi 
Ml 



(A13) 
(A14) 
(A15) 
(A16) 
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B{mlQ,Ml) 

^B{s^,MlMl) 
^B{s^,MlMl) 

B{s^,Ml,M],) 

B{s^,MlM^) 
B{s^,mlMl) 



UB{U,mlM^) 



^B{t^,mlMl) 
B{U,M^,Ml) 



B{t^,MlMl) 



A{Mf 



K) 



Ml 
Ml 



+ 



1 



+ 



1 

167r2 



1-1- 



M\ 



K 



Ind- 



Ml 



1 — (7^ In 



1-0-7, 



IGtt ' 
Ml 



+ 



1 



1 



1 A(rn|) 
2 



+ 



167r2 
1 



87r2 
1 



+ 



1 



167r2 87r2 

1 A(M|) 
Mi 



4M|- 



4M2 



\ 1/2 



arctan 



■) 



+ 



1/2 
1 

2l 



arctan 



AMI 



AMI 



(A17) 

(A18) 
(A19) 

-1/2 

1 1 (A20) 

-1/2 

(A21) 



2 ■ 167r2 

+ MkY - s^^y {nil - Mk)^ - s„ 



Ml)\n 



Ml 



X In 



IGtt^Stt 

^{me + MkY -s^ + v/(m, - Mk? 



1 A(m|) ^ 
2 



1 



2 



1 



1 



— (m 



X In 



M2 167r2 
- (m^ + M^y^t^ - {me - M^f 
^t^ - {me - M^y + ^t^ - {me + M^y 



Ml)\n(^ 



167rtvr 



^/t^ - {me - M^Y - ^t^ - {me + M^Y ' 
Vt^r - {me + M^Y^Jt^ - {me - M„)\ 



1 A{Ml) 1 A(Mi) 

2 2 M2 



167r2 



1-^(M.'-Mi)ln 



Ml 

Ml 



+ 



1 



167r2t^ 



V(M, + M;^)2 - V(M, - MkY - 



X In 



v/(M^ + M^)2 - + ^{M^ - Mk)2 - t,, 



i^r - - a/ (m^ 

2^ 1 



1 A(M2) 1 A{Ml) 

2 M^ 2 Ml 167r2 



1 - 



Mk 
1 



r2 T\/r2 



2U 



{M^ - Ml) In 



Mi 



^ ^(M, + Mk)2 - t^-^t,, - {Mr, - MkY 



X arctan 



V(M, + M^)2-i/ 
For the following integral, we shall distinguish between two cases. 



(A22) 



(A23) 



(A24) 



(A25) 



(A26) 
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(a) The lepton is an electron: 



+ 



2 Ml ' 2 M2 ' 167r2 

+\i{AMl<s^<{M,,-M^) 

v/(M, + - + v/(M, - - 
-If ((M^ - M,)2 < < {Mk - ruef) 



M2 
M2 



X In 



X arctan 



(A27) 



(b) The lepton is a muon 

B{s^,MlMl) 



l_A{M^ lA{M^ I 
2 Ml 2 167r2 



1 

2Z 



(M^ - Ml) In 



X In 



(A28) 



3. T integrals 



These integrals appeared while splitting strong and electromagnetic parts in two-point 
functions. We are interested in the following particular r integrals: 



^T{s^,MlMl) 
^T{s,,MlMl) 

r(s^, Ml, Ml) 
T{t^,M^,Ml) 



T{t^,MlMl) 



1- (T-jr 



27r 



4 fAMl 



- 1 



-1/2 



AMI 
arctan I 1 



-1/2 



v/(M^ - Mk)^ - U^y{M^ + Mk)^ - U 

X In ^ ~ ^ ~ ~ 

+ ^i^)' -t-- ViM. - Mk)^ - U 
4 

y/U - (M, - M^^{M, + Mk)^ - U 



(A29) 
(A30) 

(A31) 



(A32) 
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X arctan 



- {Mr, - MkY 



(A33) 



4. C integrals 

These are scalar three-point functions whose definition and expressions were given in the 
appendix of Ref. j^. In what follows, we sketch some of the particular cases that we need 
for the numerical evaluation of isospin breaking corrections: 



C(m^,0,m^,0,m|,M| 



^C{Mls^,Mlnl,MlMl] 



C(Mi,s^,m2,0,M|,M|) 



167r2 





( mj \ 











+ 



Ml - mj 



In 



Ml 



K)\ 5 



(A34) 
(A35) 



3277^5^(7^ 



4Li9 



+ 



1 - Q-^ 
1 + o-^ 

2 ' 

m 



47r2 



3 

41nK)-21n(^)+ln 
21n(a,) - In 



In ( ^^A^6,) 



1 + 0", 



1 1 



(A37) 



167r2 M^m^ 1 - a] 



IK 



2\n{l-alj,)--\Yi{a,K)-\ri 



\MKmJ 



In(a^x) 



^ 1 /m, 
6 2 V M 



+Li2 ( 1 
1 



me 



K 



Mk 

1 O-^A- 



+ Li2(a|^) 



Li, 1 



K 



167r2 Mi^m^ 1 - 



IK 



In 



Mi 



+ In 



mi 



Ml - m} 



■ o'eK 



(A38) 



-lnK,)--ln 



In (cTeK] 



+Li2 (^1 - 
-Lis - 
-Li2 (l - 



O'eK + Li2 i aeK 



Mk 
Mk (TeK 



me aKK 
me cTiK 



me 

T- / I 

Ll2 I i O'eKO'KK , 

me J 



Mk (Jkk 



- Li2 1 - 



me 



K 



oeKCTKK^ I , (A39) 
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•^°-.,ln.^- 



M^ + U + xi 
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{xo -Ml + mjy - XiU, mj, M^) 
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+ Xo + \^/^it^,mlMl) 



-X'/\U,mlMl)ln 



K 

Ml 



-XQ + X^I\t^,mlMl) 
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+X'/\t,,mlMl)\n^''' 
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Ml 



mj + xi + X'/\U,mj,M^ 



m 



+ xi-X^/^{U,mlMl 
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meMx 
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